X marks the spot: Accurate energies from intersecting extrapolations of continuum quantum Monte Carlo data.

We explore the application of an extrapolative method that yields very accurate total and relative energies from variational and diffusion quantum Monte Carlo (VMC and DMC) results. For a trial wave function consisting of a small configuration interaction (CI) wave function obtained from full CI quantum Monte Carlo and reoptimized in the presence of a Jastrow factor and an optional backflow transformation, we find that the VMC and DMC energies are smooth functions of the sum of the squared coefficients of the initial CI wave function and that quadratic extrapolations of the non-backflow VMC and backflow DMC energies intersect within uncertainty of the exact total energy. With adequate statistical treatment of quasi-random fluctuations, the extrapolate and intersect with polynomials of order two method is shown to yield results in agreement with benchmark-quality total and relative energies for the C2, N2, CO2, and H2O molecules, as well as for the C2 molecule in its first electronic singlet excited state, using only small CI expansion sizes.

Toward Real Chemical Accuracy on Current Quantum Hardware Through the Transcorrelated Method.

Quantum computing is emerging as a new computational paradigm with the potential to transform several research fields including quantum chemistry. However, current hardware limitations (including limited coherence times, gate infidelities, and connectivity) hamper the implementation of most quantum algorithms and call for more noise-resilient solutions. We propose an explicitly correlated Ansatz based on the transcorrelated (TC) approach to target these major roadblocks directly. This method transfers, without any approximation, correlations from the wave function directly into the Hamiltonian, thus reducing the resources needed to achieve accurate results with noisy quantum devices. We show that the TC approach allows for shallower circuits and improves the convergence toward the complete basis set limit, providing energies within chemical accuracy to experiment with smaller basis sets and, thus, fewer qubits. We demonstrate our method by computing bond lengths, dissociation energies, and vibrational frequencies close to experimental results for the hydrogen dimer and lithium hydride using two and four qubits, respectively. To demonstrate our approach's current and near-term potential, we perform hardware experiments, where our results confirm that the TC method paves the way toward accurate quantum chemistry calculations already on today's quantum hardware.

Exploiting Locality in Full Configuration Interaction Quantum Monte Carlo for Fast Excitation Generation.

In this paper, we propose an improved excitation generation algorithm for the full configuration interaction quantum Monte Carlo method, which is particularly effective in systems described by localized orbitals. The method is an extension of the precomputed heat-bath strategy of Holmes et al., with more effective sampling of double excitations and a novel approach for nonuniform sampling of single excitations. We demonstrate the effectiveness of the algorithm for a chain of 30 hydrogen atoms with atom-localized orbitals, a stack of benzene molecules, and an Fe(II)-porphyrin model complex, whereby we show an overall efficiency gain by a factor of two to four, as measured by variance reduction per wall-clock time.

xTC: An efficient treatment of three-body interactions in transcorrelated methods.

An efficient implementation for approximate inclusion of the three-body operator arising in transcorrelated methods via exclusion of explicit three-body components (xTC) is presented and tested against results in the "HEAT" benchmark set [Tajti et al., J. Chem. Phys. 121, 011599 (2004)]. Using relatively modest basis sets and computationally simple methods, total, atomization, and formation energies within near-chemical accuracy from HEAT results were obtained. The xTC ansatz reduces the nominal scaling of the three-body part of transcorrelation by two orders of magnitude to O(N5) and can readily be used with almost any quantum chemical correlation method.

Optimizing Jastrow factors for the transcorrelated method.

We investigate the optimization of flexible tailored real-space Jastrow factors for use in the transcorrelated (TC) method in combination with highly accurate quantum chemistry methods, such as initiator full configuration interaction quantum Monte Carlo (FCIQMC). Jastrow factors obtained by minimizing the variance of the TC reference energy are found to yield better, more consistent results than those obtained by minimizing the variational energy. We compute all-electron atomization energies for the challenging first-row molecules C2, CN, N2, and O2 and find that the TC method yields chemically accurate results using only the cc-pVTZ basis set, roughly matching the accuracy of non-TC calculations with the much larger cc-pV5Z basis set. We also investigate an approximation in which pure three-body excitations are neglected from the TC-FCIQMC dynamics, saving storage and computational costs, and show that it affects relative energies negligibly. Our results demonstrate that the combination of tailored real-space Jastrow factors with the multi-configurational TC-FCIQMC method provides a route to obtaining chemical accuracy using modest basis sets, obviating the need for basis-set extrapolation and composite techniques.

Transcorrelated coupled cluster methods. II. Molecular systems.

We demonstrate the accuracy of ground-state energies of the transcorrelated Hamiltonian, employing sophisticated Jastrow factors obtained from variational Monte Carlo, together with the coupled cluster and distinguishable cluster methods at the level of singles and doubles excitations. Our results show that already with the cc-pVTZ basis, the transcorrelated distinguishable cluster method gets close to the complete basis limit and near full configuration interaction quality values for relative energies of over thirty atoms and molecules. To gauge the performance in different correlation regimes, we also investigate the breaking of the nitrogen molecule with transcorrelated coupled cluster methods. Numerical evidence is presented to further justify an efficient way to incorporate the major effects coming from the three-body integrals without explicitly introducing them into the amplitude equations.

The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry.

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

TREXIO: A file format and library for quantum chemistry.

TREXIO is an open-source file format and library developed for the storage and manipulation of data produced by quantum chemistry calculations. It is designed with the goal of providing a reliable and efficient method of storing and exchanging wave function parameters and matrix elements, making it an important tool for researchers in the field of quantum chemistry. In this work, we present an overview of the TREXIO file format and library. The library consists of a front-end implemented in the C programming language and two different back-ends: a text back-end and a binary back-end utilizing the hierarchical data format version 5 library, which enables fast read and write operations. It is compatible with a variety of platforms and has interfaces for Fortran, Python, and OCaml programming languages. In addition, a suite of tools have been developed to facilitate the use of the TREXIO format and library, including converters for popular quantum chemistry codes and utilities for validating and manipulating data stored in TREXIO files. The simplicity, versatility, and ease of use of TREXIO make it a valuable resource for researchers working with quantum chemistry data.

The association between OSA and glycemic control in diabetes.

Background: Obstructive sleep apnea (OSA) is the most common sleep-realted respiratory disorder. It is frequently comorbid with cardiovascular, cerebrovascular, and metabolic diseases and is commonly observed in populations with these comorbidities. Investigators aimed to assess the effect of OSA on glycemic control in patients with diabetes. Methods: In this cross-sectional study, 266 adult patients with diabetes mellitus (DM) attending the outpatient endocrinology clinic at the Guilan University of Medical Sciences were enrolled. Patients completed a checklist that included demographic characteristics, factors, and laboratory results in addition to Berlin and STOP-BANG questionnaires to evaluate the risk of OSA. Data were analyzed by independent t-test, Mann-Whitney U test, and Chi-squared or Fisher's exact tests using the Statistical Package for the Social Sciences (SPSS) version 17. Results: A total of 266 patients with DM were enrolled in this study (34.6% males, mean age 47.00 ± 19.04 years). Based on the Berlin Questionnaire, 38.6% of all participants were at high risk of developing OSA. Based on the STOP-BANG Questionnaire (SBQ), 45.1% were at moderate and high risks. Additionally, this questionnaire showed a significant difference between low and moderate-to-severe groups regarding sex, age, body mass index (BMI), neck size, other chronic diseases, types of DM, use of insulin, Berlin Questionnaire, fasting blood sugar (FBS), and mean HbA1c. Conclusions: Based on the SBQ, our results indicated a significant relationship between OSA and glycemic control according to mean HbA1c and FBS. Therefore, by controlling the OSA, we may find a way to acheieve better glycemic control in diabetic patients.

Density Matrix Renormalization Group for Transcorrelated Hamiltonians: Ground and Excited States in Molecules.

We present the theory of a density matrix renormalization group (DMRG) algorithm which can solve for both the ground and excited states of non-Hermitian transcorrelated Hamiltonians and show applications in molecular systems. Transcorrelation (TC) accelerates the basis set convergence rate by including known physics (such as, but not limited to, the electron-electron cusp) in the Jastrow factor used for the similarity transformation. It also improves the accuracy of approximate methods such as coupled cluster singles and doubles (CCSD) as shown by recent studies. However, the non-Hermiticity of the TC Hamiltonians poses challenges for variational methods like DMRG. Imaginary-time evolution on the matrix product state (MPS) in the DMRG framework has been proposed to circumvent this problem, but this is currently limited to treating the ground state and has lower efficiency than the time-independent DMRG (TI-DMRG) due to the need to eliminate Trotter errors. In this work, we show that with minimal changes to the existing TI-DMRG algorithm, namely, replacing the original Davidson solver with the general Davidson solver to solve the non-Hermitian effective Hamiltonians at each site for a few low-lying right eigenstates, and following the rest of the original DMRG recipe, one can find the ground and excited states with improved efficiency compared to the original DMRG when extrapolating to the infinite bond dimension limit in the same basis set. An accelerated basis set convergence rate is also observed, as expected, within the TC framework.

General Analytical Nuclear Forces and Molecular Potential Energy Surface from Full Configuration Interaction Quantum Monte Carlo.

The full configuration interaction quantum Monte Carlo (FCIQMC) is a state-of-the-art stochastic electronic structure method, providing a methodology to compute FCI-level state energies of molecular systems within a quantum chemical basis. However, especially to probe dynamics at the FCIQMC level, it is necessary to devise more efficient schemes to produce nuclear forces and potential energy surfaces (PES) from FCIQMC. In this work, we derive the general formula for nuclear forces from FCIQMC, and clarify different contributions of the total force. This method to obtain FCIQMC forces eliminates previous restrictions and can be used with frozen core approximation and free selection of orbitals, making it promising for more efficient nuclear forces calculations. After some numerical checks of this procedure on the binding curve of N2 molecule, we use the FCIQMC energy and force to obtain the full-dimensional ground state PES of the water molecule via Gaussian processes regression. The new water FCIQMC PES can be used as the basis for H2O ground state nuclear dynamics, structure optimization, and rotation-vibrational spectrum calculation.

Introduction of maltodextrin nanosponges as green extraction phases: Magnetic solid phase extraction of fluoroquinolones.

In the present work, for the first time, fully green maltodextrin nanosponges (MDNSs) as extraction phases for the preparation of magnetic sorbents were introduced. All of the raw materials and synthetic methodologies used to prepare MDNSs follow the principles of green chemistry. Fluoroquinolones (FQs) (ofloxacin, ciprofloxacin, enrofloxacin, and moxifloxacin) were used as model analytes to assess the extraction performance of MDNSs in magnetic solid-phase extraction (MSPE) procedure. Under the optimized conditions, the developed MSPE was coupled with HPLC-UV to sensitively quantify trace levels of four FQs in tap water, wastewater, and river water samples. The proposed MSPE-HPLC-UV method provided wide linearity in the range of 0.25-500.0 ng mL-1 with determination coefficients ≥0.9997. Low limits of detection (LODs ≤0.06 ng mL-1) and quantification (LOQs ≤0.20 ng mL-1) were achieved. Furthermore, the developed green extraction phases compared with all other maltodextrin-based extraction phases were promising.

Importance-sampling FCIQMC: Solving weak sign-problem systems.

We investigate the exact full configuration interaction quantum Monte Carlo algorithm (without the initiator approximation) applied to weak sign-problem fermionic systems, namely, systems in which the energy gap to the corresponding sign-free or "stoquastized" state is small. We show that the minimum number of walkers required to exactly overcome the sign problem can be significantly reduced via an importance-sampling similarity transformation even though the similarity-transformed Hamiltonian has the same stoquastic gap as the untransformed one. Furthermore, we show that in the off-half-filling Hubbard model at U/t = 8, the real-space (site) representation has a much weaker sign problem compared to the momentum space representation. By applying importance sampling using a Gutzwiller-like guiding wavefunction, we are able to substantially reduce the minimum number of walkers in the case of 2 × â„“ Hubbard ladders, enabling us to get exact energies for sizable ladders. With these results, we calculate the fundamental charge gap ΔEfund = E(N + 1) + E(N - 1) - 2E(N) for the ladder systems compared to strictly one-dimensional Hubbard chains and show that the ladder systems have a reduced fundamental gap compared to the 1D chains.

Perturbation calculation of the uniform electron gas with a transcorrelated Hamiltonian.

With a transcorrelated Hamiltonian, we perform a many body perturbation calculation on the uniform electron gas in the high density regime. By using a correlation factor optimized for a single determinant Jastrow ansatz, the second order correlation energy is calculated as 1-ln⁡2π2ln(rs)-0.05075. This already reproduces the exact logarithmic term of the random phase approximation (RPA) result, while the constant term is roughly 7% larger than the RPA one. The close agreement with the RPA method demonstrates that the transcorrelated method offers a viable and potentially efficient method for treating metallic systems.

Performance of a one-parameter correlation factor for transcorrelation: Study on a series of second row atomic and molecular systems.

In this work, we investigate the performance of a recently proposed transcorrelated (TC) approach based on a single-parameter correlation factor [E. Giner, J. Chem. Phys. 154, 084119 (2021)] for systems involving more than two electrons. The benefit of such an approach relies on its simplicity as efficient numerical-analytical schemes can be set up to compute the two- and three-body integrals occurring in the effective TC Hamiltonian. To obtain accurate ground state energies within a given basis set, the present TC scheme is coupled to the recently proposed TC-full configuration interaction quantum Monte Carlo method [Cohen et al., J. Chem. Phys. 151, 061101 (2019)]. We report ground state total energies on the Li-Ne series, together with their first cations, computed with increasingly large basis sets and compare to more elaborate correlation factors involving electron-electron-nucleus coordinates. Numerical results on the Li-Ne ionization potentials show that the use of the single-parameter correlation factor brings on average only a slightly lower accuracy (1.2 mH) in a triple-zeta quality basis set with respect to a more sophisticated correlation factor. However, already using a quadruple-zeta quality basis set yields results within chemical accuracy to complete basis set limit results when using this novel single-parameter correlation factor. Calculations on the H2O, CH2, and FH molecules show that a similar precision can be obtained within a triple-zeta quality basis set for the atomization energies of molecular systems.

FCIQMC-Tailored Distinguishable Cluster Approach: Open-Shell Systems.

A recently proposed tailored approach based on the distinguishable cluster method and the stochastic FCI solver, FCIQMC [J. Chem. Theory Comput. 2020, 16, 5621], is extended to open-shell molecular systems. The method is employed to calculate spin gaps of various Fe(II) complexes, including a Fe(II) porphyrin model system. Both distinguishable cluster and fully relaxed CASSCF natural orbitals were used in this work as reference for the subsequent tailored distinguishable cluster calculations. The distinguishable cluster natural orbitals occupation numbers were also used as an aid to the selection of the active space. The effect of the active space sizes and of the explicit correlation correction (F12) onto the predicted spin gaps is investigated. The tailored distinguishable cluster with singles and doubles yields consistently more accurate results compared to the tailored coupled cluster with singles and doubles.

Association Between Poor Sleep Quality and Glycemic Control in Adult Patients with Diabetes Referred to Endocrinology Clinic of Guilan: A Cross-sectional Study.

Background: Diabetes is a prevalent chronic medical comorbid condition worldwide. Diabetes mellitus is associated with various sleep disorders. Objectives: We aimed to determine the prevalence of poor sleep and the main factors of sleep interruptions in patients with diabetes mellitus. We further evaluated the association of sleep interruptions with glycemic control in this cohort. Methods: We conducted a cross-sectional study on 266 patients with type 1 and type 2 diabetes recruited from a university outpatient endocrinology clinic. Patients completed a checklist including demographic and disease-related characteristics in addition to the Pittsburgh Sleep Quality Index (PSQI) to evaluate sleep quality. Using the PSQI cutoff score of 5, we created two subgroups of good sleepers (GS) and poor sleepers (PS). Results: Our results showed that good sleeper and poor sleeper patients with diabetes were significantly different regarding sex, employment status, BMI, presence of diabetes-related complications, HbA1c, and 2-hour postprandial blood sugar (2HPPBS) (all significant at P < 0.05). The most prevalent factors of sleep interruptions were "waking up to use a bathroom", "feeling hot", "pain", "having coughs or snores", and "bad dreams". Among the subjective factors of sleep interruption, problems with sleep initiation, maintenance, or early morning awakenings in addition to having pain or respiratory problems such as coughing or snoring had the most significant associations with HbA1c. Conclusions: Our study showed significant subjective sleep disturbances (both quality and quantity) in patients with diabetes mellitus (both type I and II) and its association with diabetes control. We further identified the main factors that led to sleep interruptions in this cohort.

Full configuration interaction quantum Monte Carlo treatment of fragments embedded in a periodic mean field.

We present an embedded fragment approach for high-level quantum chemical calculations on local features in periodic systems. The fragment is defined as a set of localized orbitals (occupied and virtual) corresponding to a converged periodic Hartree-Fock solution. These orbitals serve as the basis for the in-fragment post-Hartree-Fock treatment. The embedding field for the fragment, consisting of the Coulomb and exchange potential from the rest of the crystal, is included in the fragment's one-electron Hamiltonian. As an application of the embedded fragment approach, we investigate the performance of full configuration interaction quantum Monte Carlo (FCIQMC) with the adaptive shift. As the orbital choice, we use the natural orbitals from the distinguishable cluster method with singles and doubles. FCIQMC is a stochastic approximation to the full CI method and can be routinely applied to much larger active spaces than the latter. This makes this method especially attractive in the context of open shell defects in crystals, where fragments of adequate size can be rather large. As a test case, we consider dissociation of a fluorine atom from a fluorographane surface. This process poses a challenge for high-level electronic structure models as both the static and dynamic correlations are essential here. Furthermore, the active space for an adequate fragment (32 electrons in 173 orbitals) is already quite large even for FCIQMC. Despite this, FCIQMC delivers accurate dissociation and total energies.

General embedded cluster protocol for accurate modeling of oxygen vacancies in metal-oxides.

The O vacancy (Ov) formation energy, EOv, is an important property of a metal-oxide, governing its performance in applications such as fuel cells or heterogeneous catalysis. These defects are routinely studied with density functional theory (DFT). However, it is well-recognized that standard DFT formulations (e.g., the generalized gradient approximation) are insufficient for modeling the Ov, requiring higher levels of theory. The embedded cluster method offers a promising approach to compute EOv accurately, giving access to all electronic structure methods. Central to this approach is the construction of quantum(-mechanically treated) clusters placed within suitable embedding environments. Unfortunately, current approaches to constructing the quantum clusters either require large system sizes, preventing application of high-level methods, or require significant manual input, preventing investigations of multiple systems simultaneously. In this work, we present a systematic and general quantum cluster design protocol that can determine small converged quantum clusters for studying the Ov in metal-oxides with accurate methods, such as local coupled cluster with single, double, and perturbative triple excitations. We apply this protocol to study the Ov in the bulk and surface planes of rutile TiO2 and rock salt MgO, producing the first accurate and well-converged determinations of EOv with this method. These reference values are used to benchmark exchange-correlation functionals in DFT, and we find that all the studied functionals underestimate EOv, with the average error decreasing along the rungs of Jacob's ladder. This protocol is automatable for high-throughput calculations and can be generalized to study other point defects or adsorbates.

Spin Purification in Full-CI Quantum Monte Carlo via a First-Order Penalty Approach.

In this article, we demonstrate that a first-order spin penalty scheme can be efficiently applied to the Slater determinant based Full-CI Quantum Monte Carlo (FCIQMC) algorithm, as a practical route toward spin purification. Two crucial applications are presented to demonstrate the validity and robustness of this scheme: the 1Δg ← 3Σg vertical excitation in O2 and key spin gaps in a [Mn3(IV)O4] cluster. In the absence of a robust spin adaptation/purification technique, both applications would be unattainable by Slater determinant based ground state methods, with any starting wave function collapsing into the higher-spin ground state during the optimization. This strategy can be coupled to other algorithms that use the Slater determinant based FCIQMC algorithm as configuration interaction eigensolver, including the Stochastic Generalized Active Space, the similarity-transformed FCIQMC, the tailored-CC, and second-order perturbation theory approaches. Moreover, in contrast to the GUGA-FCIQMC technique, this strategy features both spin projection and total spin adaptation, making it appealing when solving anisotropic Hamiltonians. It also provides spin-resolved reduced density matrices, important for the investigation of spin-dependent properties in polynuclear transition metal clusters, such as the hyperfine-coupling constants.